D’THfjEY TITr'T T r nv ’ ». HY 

F n » P^HOOL 



NON -NEWTONIAN TEMPERATURE AND PRESSURE EFFECTS OF A POWDER 
LUBRICANT SLURRY IN A ROTATING HYDROSTATIC STEP BEARING 



BY 

WILLIAM EDWARD FINN 

/ 



MASTER OF SCIENCE IN 
MECHANICAL ENGINEERING 



A THESIS PRESENTED TO THE GRADUATE SCHOOL 
OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT 
OF THE REQUIREMENTS FOR THE DEGREE OF 
MASTER OF SCIENCE 

UNIVERSITY OF FLORIDA 

1993 



T 25971 3 



DTXDLE^ K?VT- r .TT;r: ,*,T7 
HAVAL PG3TG t. * V new nr 
MOKTEAE ¥. Cit-EiEGHEiA 9<aHvo out/*' 



ACKNOWLEDGMENTS 



I am very grateful for the support provided by the 
following individuals: 

Dr. D.W. Dareing first interested me in the tribology 
area and was instrumental in guiding my research. He helped 
refine the initial problem statement that became my master's 
thesis . 

Dr. J.E. Peterson devoted many hours obtaining funding, 
providing counsel, solving experimental setbacks and 
supervising my research. Her encouragement and assistance 
were critical to the timely completion of this thesis. 

Dr. J.K. Schueller's motivation and interest kept the 
project moving and provided additional reference material 
for this project. 

Finally, the test rig modifications would not have been 
possible without the technical assistance from Mr. Steve 
Sowa and his division's machining work, material from Mr. 
Howard Purdy and Mr. Jeff Studstill and electrical guidance 
and material from Mr. Robert Harker. 



ii 



TABLE OF CONTENTS 



ACKNOWLEDGMENTS ii 

NOMENCLATURE v 

ABSTRACT vii 

CHAPTERS 

1 BACKGROUND AND SIGNIFICANCE OF RESEARCH. . . 1 

Background of Hydrostatic Bearings 1 

Significance of the Research 2 

2 OBJECTIVE OF THE RESEARCH 4 

3 RESEARCH EQUIPMENT AND METHOD 6 

Laboratory Test Equipment 6 

Research Method 9 

Uncertainty Analysis and Calibration 10 

4 MATHEMATICAL MODELING 17 

Temperature' Equations 19 

Ethylene Glycol Relations 19 

Graphite/Ethylene Glycol Relations 22 

Pressure Equations 27 

Ethylene Glycol Relations 28 

Graphite/Ethylene Glycol Relations 30 



iii 



5 EXPERIMENTAL RESULTS 32 

Experimental Method 32 

Temperature Results 32 

Pressure Results 35 

6 CONCLUSIONS 42 

REFERENCES 45 

APPENDIX SAMPLE UNCERTAINTY CALCULATION 4 7 

BIOGRAPHICAL SKETCH 4 8 



NOMENCLATURE 



A Area of bearing surface, in 2 

A, B... Linear regression coefficients 
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V Kinematic fluid viscosity, in 2 /sec 

co Angular velocity, rad/sec 
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The purpose of this research was to investigate the 
pressure and temperature effects of graphite powder 
lubricant added to a carrier fluid of ethylene glycol as 
applied in a rotating hydrostatic step bearing. After the 
testing apparatus was designed and modified, the temperature 
and pressure profiles were determined analytically and 
experimentally . 

The rheological behavior of the non -Newtonian lubricant 
was modeled using a Power Law model developed previously. 
This model has been proven to match experimental data and 
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was therefore the tool used for comparison with research 
data. 

Ethylene glycol was used as the Newtonian lubricant for 
the initial experimentation, which provided a check on the 
test apparatus and a comparison with the non-Newtonian 
graphite slurry. 

The measured data for the temperature increase 
dependency on rotation for both fluids compared favorably to 
the mathematical predictions. A significantly higher 
temperature difference was seen in the non-Newtonian 
lubricant due to the higher shear rates. 

The pressure profile, while not directly dependent on 
rotation as shown in the mathematical model, demonstrated a 
reduction based upon higher rotation speeds. This reduction 
was greater for the non-Newtonian lubricant. Due to the 
higher temperatures seen in the non-Newtonian lubricant and 
temperature dependent Power Law constants, a pressure loss 
is expected. For the Newtonian lubricant, viscosity is 
dependent upon fluid temperature. 

Therefore, the effects of operating speed and 
temperature on a non-Newtonian lubricant should be 
considered, as well as their greater load-carrying capacity. 
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CHAPTER 1 

BACKGROUND AND SIGNIFICANCE OF RESEARCH 

Background of Hydrostatic Bearings 
Hydrostatic bearings were first demonstrated in 1851 by 
L.D. Girard through an invention involving high pressure 
water- fed bearings used for a railway propulsion system and 
patented in 1865 (1, p. 225; 2, p. 435). This externally 
pressurized or hydrostatic system was further demonstrated 
at the 1878 Paris Industrial Exposition when a heavy block 
with four feet was floated on a thin film of lubricant that 
was pumped down each leg (1, 2) . The first analysis of a 
thrust bearing was provided by Lord Rayleigh in 1917 using 
the hydrostatic system and solving the load, flowrate and 
frictional torque (1; 2, p. 436) . 

Current applications include the analysing magnet of 
the Nuclear Structure Facility at the SERC Daresbury 
Laboratory (l, p. 227) , the Halle optical telescope (2, p. 
436), the 210 - foot-diameter tracking antenna at the NASA-Jet 
Propulsion Laboratory Deep Space Instrumentation Facility 
(3, p. 33-51), and Denver's Mile High Stadium (4, p. 78) . 
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All applications involve moving extreme weight with minimal 
effort and in some cases, high precision. In 1972, the 
Institution of Mechanical Engineers published proceedings on 
a conference solely devoted to externally pressurized 
bearings (2, p. 437) . Since that time, more effort has been 
spent by industry on the study of the lubricants as most of 
the basic mechanical principles of the hydrostatic bearing 
had been solved. 

Significance of the Research 
The addition of solid lubricants to Newtonian fluids 
has been the subject of extensive research. One major 
reason for the addition of solid lubricants is the greater 
stability of these mixtures when subjected to extreme 
temperatures. The aerospace industry, military, and NASA 
all report using solid lubricants because of their greater 
load and temperature capabilities (5, p. 5; 6, p. 5, 83-86) . 
Many commercial applications exist for a Graphite/Glycol 
lubricant mixture including switch gear equipment and high 
temperature oven bearings (5, p. 85-86). 

The U.S. gas turbine industry has started a new program 
in cooperation with NASA and the Department of Defense 
called the Integrated High Performance Turbine Engine 
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Technology (IHPTET) program (7, p. 114) . The desired double 
propulsion capability from the large speeds in the gas 
turbine shafts result in high temperatures which have 
further emphasized the requirement for nontraditional 
lubricants. The approximate working temperatures of 
different types of lubricants are as follows (8, p. 2) : 
Petroleum Lubricants below 100°C 

Synthetic Lubricants 200 to 315°C 

Powder Lubricants 425 to 675°C 

The powder lubricant is supplied to the system through 
pressurized air or slurries containing the powder (7, p. 

114) . These non-Newtonian lubricants are the subject of 
current research towards more accurate behavior prediction. 



CHAPTER 2 

OBJECTIVE OF THE RESEARCH 



This research investigates the pressure and temperature 
affects of graphite powder lubricant added to a carrier 
fluid of ethylene glycol as applied in a rotating 
hydrostatic step bearing. The graphite powder was provided 
by Union Carbide and the ethylene glycol was provided by 
Fisher Scientific. The methodology included laboratory 
testing and mathematical modeling for comparison with the 
experimental results. 

An existing laboratory test rig was extensively 
modified to impart an angular velocity to a hydrostatic step 
bearing and to measure temperature, pressure, flowrate, 
revolutions per minute (RPM) , and film thickness. 
Measurements for the ethylene glycol were obtained for use 
as a check on the testing rig and for establishing baseline 
data . 

The goal was to formulate temperature equations 
dependent upon the rotational velocity and compare this 
affect with the pressure distribution across the bearing 
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radius. Potential inertial affects from the step bearing's 
rotation and its effects on the graphite particles in the 
ethylene glycol would also be examined in the pressure 
profile . 

An understanding of these fundamental lubrication 
principles in a non-Newtonian fluid is necessary for 
developing a powder lubricant slurry approach in future 



applications . 



CHAPTER 3 

RESEARCH EQUIPMENT AND METHOD 



Laboratory Test Equipment 

The laboratory testing was an important portion of this 
research. Significant time and effort was spent designing, 
fabricating, and redesigning the modifications to the 
existing test rig. The modified hydrostatic test rig 
(Figure 3-1) contains a flat base and a parallel hydrostatic 
step bearing (items 5 and 6, Figure 3-1) . The bearing 
consisted of an upper surface with an outside diameter of 
four inches and an inside diameter of two inches. It was 
connected to the shaft and gear assembly through a ball 
bearing assembly designed and fabricated to maintain a 
parallel surface to the base while rotating. The resultant 
constant film thickness was measured by three LVDT 
displacement sensors (item 4, Figure 3-1) mounted on 
supports which were placed 120 degrees apart on top of the 
hydrostatic step bearing. Since three points define a 
plane, the use of three sensors was sufficient. The voltage 
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Figure 3-1. Modified Hydrostatic Test Rig 
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reading of the LVDT displacement sensors was compared to the 



calibration and provided film thickness. The specifications 



for the LVDT sensors are as follows: 



Make 

Model 

Input Volts DC (v) 
Linearity % 

Freq. Response 3dB down 
Temperature Range 
Resolution 



Trans -Tek, Inc. 
DC-DC Series 240 
02141-0000 
6 to 30 
.5 

140 Hz 

-50 C to 120 C 
Infinite 



The film thickness could be adjusted through the fine 



adjusting thread (item 1, Figure 3-1) that connected the 



screw shaft to the frame (item 2, Figure 3-1) . The gear 



assembly was belt driven by a controllable, variable speed 
motor (item 3, Figure 3-1) . The RPM was measured directly 



through a Pioneer Digital Photo tachometer, model DT-36. 



The motor specifications are as follows: 



Make 

Model 

HP 

RPM 

Armature Volts 
Armature Current 



Dayton Permanent 
Magnet DC 
2M168B 
1/2 
1725 
90 VDC 
5 . 5 amps 



Two thermisters were installed in the base of the 



bearing surface (item 7, Figure 3-1) and measured the 
temperature of the lubricant entering the recess pocket and 
at the perimeter of the hydrostatic step bearing. The 
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measured resistance was compared to the thermister's 
calibration and provided the temperature. 

The lubricant was supplied to the bearing surfaces by 
an external pump through an internal orifice. The pump's 
flow rate was adjustable. It's specifications are as 



follows : 



Make 

Model 

Mo tor /Pump Coupling 
Maximum Flow Rate 
Rated Maximum Pressure 



Micropump Corp. 
415 

Magnetic Drive 
630 ml /min 
20 psi 



The fluid flow was routed through a rotameter flow 



meter. Comparing the reading with the flow meter's 



calibration provided the flow rate. 



The fluid film pressures were measured with dial 
pressure gages at four locations across the bearing surface 
located: 



1) within the recess pocket 

2) 1/4 inch from the inside pocket 

3) 1/2 inch from the inside pocket 

4) 3/4 inch from the inside pocket 

Research Method 



The experimental test rig was used to determine the 
effects of powder lubricants in a carrier fluid subjected to 
rotation in a hydrostatic step bearing. The testing method 
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is described in Chapter 5. First, the various measuring 
devices were calibrated. 

Uncertainty Analysis and Calibration 
The uncertainty analysis for the flowmeter, thermister, 
and LVDT data was calculated using a least squares fit 
technique to determine the deviation width. Since each 
curve fit was linear, the formula (9, p. 158) used to 
determine the standard deviation is written as equation 
(3-1) where the coefficients A and B are determined through 
the linear regression curve fit. Further, <j y was multiplied 
by 1.7 in order to obtain 10:1 odds or a 90% probability 
that a measurement will fall within 1.7 standard deviations. 

vj = -^'L(y i -A-Bxi) 2 (3-D 

The rotameter flowmeter was calibrated by adjusting the 
flowrate through the pump, marking the flowmeter reading, 
and timing the fill time of the liquid into a beaker of 
known volume. A linear regression or least- squares fit (9, 
p. 154) for a line of the measured data provided the 



calibration equation: 
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y = 5 . 55 * - 205.69 



(3-2) 



where y is the flow rate in ml/min and x is the 
flowmeter reading. The calibration curve is shown in figure 
3-2. The flowmeter uncertainty calculated through equation 
(3-1) for both liquids was +11.8 ml/min. 




* Data Eq. (3-2) 



Figure 3-2. Flowmeter calibration curve 



The thermisters were calibrated by placing them in a 
water bath, adjusting the temperature, reading the mercury 
thermometer temperature, and recording the thermisters' 
resistance. The uncertainty of the thermometer using half 
the least count was ±0.1°F. Thermisters function on a 
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principle that their resistance is lower at higher 
temperatures. This dependency is logarithmic and shown by 
the equation: 



7 ’ 0 ) + ^ {3 ' 3) 

where R 0 and T 0 are the baseline resistance in kOhms 

and temperature in Rankine, respectively. The constants A 
and B were determined through calibration for the two 
thermisters using linear regression. Rearranging the 
results and solving for resistance as a function of 
temperature provided the calibration equations: 



R i = 11.247 

f 

R 2 = 11.56 



' io 3083 - 39 fe -^)-° 0092 

10 3ll2 ' 86 fe-^h 00027 ^ 



\ 

J 






J 



(3-4) 

(3-5) 



where R x , T x and R 2/ T 2 are the resistances and 
temperatures of thermisters one and two, respectively. The 
calibration curves are shown in figures 3-3 and 3-4. Their 
uncertainty calculated through equation (3-1) was ± 1 Ohm 



( . 001 kOhms) . 
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* Thermister #1 Eq. (3-4) 



Figure 3-3. Thermister One calibration curve 




° Thermister #2 Eq. (3-5) 



Figure 3-4. Thermister Two calibration curve 
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The LVDTs were calibrated by unscrewing the shaft of 
known adjusting thread through equidistant intervals 
beginning at zero when the hydrostatic step bearing was 
touching the flat base. The voltage reading of each LVDT 
was recorded at these intervals. Linear regression of the 
measured data for each LVDT provided the following 
calibration equations: 

LVDT A: = -1 4S.4Xa + 1 9.44 (3-6) 

LVDT B: +# = —1 33.8*8 + 1 9.78 (3-7) 

lvdt c: yc = ~97. 6xc+ 19.97 ( 3 - 8 ) 

where y is the LVDT reading in volts and x is the film 
thickness in inches. The calibration curves are shown in 
figure 3-5. The uncertainty calculated through equation 
(3-1) was ±0.04 volts for LVDT A, +0.07 volts for LVDT B, 
and +0.05 volts for LVDT C. 

The uncertainty for the tachometer was ±1 RPM as given 
by the manufacturer. It was verified by repeatedly setting 
the maximum RPM value and noting that the value was 
consistently within ±1 of 670 RPM. The uncertainty for the 
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Film Thickness (inches) 



* LVDT A 


U LVDT B 


* LVDT C 


Eq. (3-6) - 


Eq. (3-7) 


Eq. (3-8) 



Figure 3-5. LVDT calibration curves 



dial pressure gages was +0.5 psi using half the least count. 

The film thickness, h, was calculated as an average of 
the three LVDT readings . In order to determine the overall 
uncertainty of the film thickness during the testing, 
equations (3-6), (3-7), and (3-8) were rearranged to solve 

for film thickness based upon LVDT voltage. The three 
equations were averaged resulting in 

h = = -.0022x4 - ,0025y B - .0033>>c + . 1 6 1 < 3 - 9 > 
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A root- sum- square method (9, p. 73; 10, p. 15; 11, p. 
254) was employed as shown in equation (3-10) where h is 



determined from equation (3-9) 



. dh 
and ^r“ 



is its partial 



derivative with respect to LVDT A, B, or C. The uncertainty 
of LVDT A, B, or C is COi as derived earlier. 



0>A = 



(«-)■' 



1/2 



(3-10) 



The overall uncertainty of h, the film thickness, was 
±0.0003 inches. This is important because the film 
thickness in the mathematical modeling of the pressure 
profile is raised to the third power resulting in a three 
times greater uncertainty. 



CHAPTER 4 

MATHEMATICAL MODELING 



The mathematical calculations in this chapter were 
formulated in order to model the temperature change and the 
pressure profile in a hydrostatic step bearing while varying 
the film thickness and revolutions per minute (RPM) . The 
comparisons of the measured data with the mathematical 
modeling are discussed in Chapter 5. 
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The equations derived for the temperature and pressure 
relationships of the hydrostatic step bearing shown in 
figure 4-1 involve some basic assumptions (4, pp. 80-88; 8, 

p. 10) : 

1. The flow is laminar. The maximum Reynolds number 
using the greatest flow rate and the radius at the 
step is 25.9. A Reynolds number less than 1000 in 
bearing applications is assumed to be laminar flow 
(4, p. 63; 12, p. 343-346) . The flow is pressure 
induced . 

2. The inertial forces resulting from the rotation of 
the bearing are negligible when compared to the 
viscous shear forces and can be neglected. 

3. The fluids are incompressible. 

4. The pressure is a function of the radius, 
neglecting any side leakage. 

5. The temperature is constant across the film 
thickness . 

6. All the heat generated appears in the fluid and 
corresponds to its temperature rise. 

7. The pressure at the outside of the bearing is 



assumed to be ambient. 
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Temperature Equations 

Because of the different properties of Newtonian and 
non-Newtonian fluids, two different mathematical model 
relations were used. 

Ethylene Glycol Relations 

Pure ethylene glycol is a Newtonian fluid and 
therefore, the following relations apply according to the 
geometry of the test rig as shown in the schematic (figure 
4-1) . 

The pump loss is given by the relation (4, p. 86), 

E p = QAP (4-i) 

where Q is the flow rate in cubic inches per second, 

Ap is the pressure difference in psi, and E p is energy in 
inch-pounds per second. The friction loss on the bearing 
due to rotation can be solved beginning with Newton's fluid 
friction law, 






(4-2) 
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Referring to Figure 4-1 and substituting the film 
thickness, h, for the differential width, this is equivalent 
to Newton's equation of viscous friction drag. 






(4-3) 



Substituting nr 2 for A and considering a differential 
annulus results, 

dF = irdr) ( 4 - 4 ) 



When one of the contacting surfaces is stationary, 
v=cor, where CD is the angular velocity in radians per 
second. Multiplying both sides by r to determine the torque 
on an elemental annulus of radius r and width dr is 



rdF = dM— ^j^-r 3 dr 

h 



The total frictional torque (4, p. 88) becomes 



M = r 3 dr = 

h J R„ 



R 

Ro 




(4-6) 
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where |i is the viscosity in reyns, to is the angular 
velocity in radians per second, h is the film thickness in 
inches, R 0 is the radius to the step in inches, and R is the 
radius to the outside edge of the bearing in inches. 

The friction loss is given by the equation (4, p. 88) 




2 nMN 
60 



(4-7) 



where E f is energy in inch-pounds per second, N is the 
rotation in revolutions per minute (RPM) , and M is the 
torque in inch-pounds. Therefore, the total energy loss 
across the bearing from r=o to r=R can be expressed as 
follows (4, p. 90) , 



Et = E p +Ef= \2JQyCAT (4-8) 

where 12J is the conversion factor for foot-pounds to 
British thermal units (Btu) , Q is the flow rate of the fluid 
in cubic inches per second, y is the weight density of the 
fluid in pounds per cubic inches, C is the specific heat of 
the fluid in Btu per pounds- °F, and AT is the temperature 
difference across the bearing in °F. The temperature 
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difference can be solved for as a function of the rotation 
since the flow rate, pressure difference, and film thickness 
are held constant. The pressure difference from r=0 to r=R 
equals the supply pressure since the pressure at the outside 
of the bearing, P(R), is assumed to be ambient. Therefore, 
substituting the values for E f and E p and rearranging to 
solve for AT results 



Graphite /Ethylene Glycol Relations 

A mixture of ethylene glycol and graphite powder is a 
non-Newtonian fluid and requires different fundamental 
modeling. Currently, there are three accepted rheological 
models (7, pp. 115-116; 8, p. 15) . First, the Bingham model 
is stated: 




(4-9) 



Non-N — + K 0 y 



(4-10) 



where x Non . N is the non-Newtonian shear stress in 



Pascals, T 0 is the yield point in Pascals, K 0 is the Bingham 
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constant in Pascal - seconds , and y is the shear rate in 
seconds' 1 . This model has shown to deviate significantly 
from measurements at low shear rates, however, it is 
accurate at high shear rates where many lubricant films 
operate (7, p. 115) . Second, the Hyperbolic model is 
stated: 



where the shear stress, yield point, Bingham constant, 
and shear rate variables are the same as the Bingham model. 
This function is similar to the Bingham model at high shear 
rates and has proven to match better with results at low 
shear rates (7, p. 116) . 

Finally, the Power Law model is stated: 



where x Non . N is the non-Newtonian shear stress in 
Pascals, K is the Power Law constant in Pascal - seconds", y 
is the shear rate in seconds", and n is a nondimens ional 
Power Law constant. This model has proven useful for 




(4-11) 



^ Non-N — Ky n 



(4-12) 
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describing shear stress -shear rate relations for 
non-Newtonian fluids without a yield point and is useful for 
fluid flow formulations (7, p. 116) . 

Since it has been shown that the Power Law model can 
reasonably model a 1:8 ratio of graphite powder to ethylene 
glycol (7, p. 116; 8, pp. 31-32; 13, pp. 86-87) and the 
mathematical expression is uncomplicated to evaluate, it was 
used in this research. 



Table 4-1. Power Law Constants 



Graphite -Ethylene Glycol (1:8) 


Temp, °C 


K, Pa- sec" 


n 


20 


0 . 0572 


0.973 


30 


0 . 0448 


0.971 


40 


0.039 


0 .969 


50 


0.0361 


0.969 



The Power Law constants used were taken from Table 4-1 
(7, p. 116) dependent on the temperature of the fluid. 

The pump loss given by equation (4-1) is still 
applicable. The friction loss on the bearing due to 
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rotation can be determined beginning with the Power Law 
model, equation (4-12), rewritten as follows: 

'= K {%)" (4 - 13 > 

This is equivalent to the equation of viscous friction 

drag, 



F = KA 




(4-14) 



Substituting jtr 2 for A and considering a differential 
annulus results, 



dF = K{l 



y 

J 



(2nrdr) 



(4-15) 



Since one of the contacting surfaces is stationary, 
v=cor, where co is the angular velocity in radians per 
second. Multiplying both sides by r to determine the torque 
on an elemental annulus of radius r and width dr results, 



rdF = dM= 2%K^jj (r n+1 dr) 



(4-16) 
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The total frictional torque becomes 




(4-17) 



This is equivalent to equation (4-6) for the Newtonian 
fluid where n=l and K is the viscosity. A conversion factor 
is necessary since the Power Law constant, K, has units of 
Pascals versus pounds per square inch (psi) . The friction 
loss from equation (4-7) expands to the following: 



Therefore, the total energy loss from r=0 to r=R can be 
determined according to equation (4-8) . The temperature 
difference is solved as a function of the rotation since the 
flow rate, pressure difference, and film thickness are held 
constant. The pressure difference from r=0 to r=R equals 
the supply pressure since the pressure at the outside of the 




(4-18) 



bearing, P(R), is assumed to be ambient. Therefore, 
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substituting the values for E f and E p and rearranging to 
solve for At results 



AT= 



n 2 KN 



12/yC 1240200 /£>yC^30/! 




R n+2 —Ko + ^ 
n + 3 



] 



(4-19) 



Pressure Equations 

The different properties for Newtonian and 
non-Newtonian fluids required different mathematical model 
relations for each as previously discussed in the 




Figure 4-2. Pressure distribution on a unit volume 
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derivations of equations (4-2) and (4-12) . The fluid 
flow in the hydrostatic step bearing is pressure induced. 
Ethylene Glycol Relations 

For the geometry of the hydrostatic step bearing 
(Figure 4-1) , the pressure distribution on a unit volume is 
shown in Figure 4-2. 

Summing the forces in the horizontal direction and 
rearranging is 



dP_dx 
dx dy 



(4-20) 



Substituting equation (4-2) , it is a simplified form of 
the Navier- Stokes equations (12, p.343) for pressure induced 
flow between two parallel surfaces. 



— = | | d 2y 

dx ~ 



(4-21) 



Integrating twice across the y direction, assuming 
dp/dx is a constant across the film thickness, yields 

v M = i:^y 2+c iy +c 2 



(4-22) 
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Evaluating the constants of integration, C 1 and C 2 , 
with the boundary conditions that v(0)=0 and v(h)=0 results 



Integrating, it is stated (4, p. 81; 7, p. 118; 8, p. 

12 ) 



where Q is the flow rate in cubic inches per second, b 
is the plate width in inches, h is the film thickness in 
inches, |i is the viscosity in reyns, and dp/cbc is the 
pressure gradient. Equation (4-25) is commonly used in 
radial flow situations. Since ethylene glycol is a 
Newtonian fluid, Replacing b by 2;tr, dx by dr, integrating 
and solving for the pressure at any radius, r, becomes (4, 




(4-23) 



Solving for the flow rate, Q, is shown 



Q = \ o v{y)dyb 



(4-24) 



q _ bP dp 
12(1 etc 



(4-25) 



p. 82; 8, p . 14 ; 14, p. 7) 
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(4-26) 



where P is the pressure at any point along the radius 
in pounds per square inch (psi) , R is the outside radius of 
the bearing, r is any radius between R 0 and R, Q is the flow 
rate in cubic inches per second, h is the film thickness in 
inches, and ji is the viscosity in reyns. 

Graphite /Ethylene Glycol Relations 

Using the same assumptions as the ethylene glycol 
relations with the exception of equation (4-12) for the 
Power Law Model described earlier, it is stated (7, p. 118; 
8, p. 16; 14, p. 7) 



Replacing b by 2izr, dx by dr, integrating and solving 
for the pressure at any radius, r, becomes (8, p. 16; 14, p. 
8 ) 




(4-27) 




(4-28) 
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The Power Law constants, K and n, are determined from 
table 4-1 for equations (4-22) and (4-23) dependent upon the 
temperature of the fluid. 

For this research, the pressure values for both fluids 
were nondimensionalized by dividing by the inlet supply 
pressure. The flow rate and film thickness were constant so 
that the RPM could be adjusted in order to observe the 
rotational effects on the pressure at the set points along 
the radius . 



CHAPTER 5 

EXPERIMENTAL RESULTS 



Experimental Method 

The initial testing with the laboratory test rig was 
performed with pure ethylene glycol as the lubricant. This 
data provided a check on the instruments and equipment and a 
comparison for studying the effects of adding graphite 
powder to the ethylene glycol. 

The film thickness was set based upon the LVDT 
readings. The flow rate of the lubricant was adjusted so 
that the maximum pressure of the dial gages would not be 
exceeded. The motor was energized and readings were taken of 
the LVDT voltages, flowmeter, pressure gages, and 
thermisters as the rotation was adjusted from 100 to 670 
RPM. When the Graphite/Ethylene Glycol slurry was the 
lubricant, equivalent film thicknesses as the pure ethylene 
glycol were obtained and the above procedure was performed. 

Temperature Results 

The purpose of the thermister data was to determine the 
relation between the temperature difference at the inlet and 
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the edge of the bearing and the rotation speed (RPM) . This 
experimental data were then compared with the mathematical 
predictions as developed in Chapter 4. The data shown in 
Figures 5-1 and 5-2 closely match the predictions within the 
uncertainty. The repeatability was within experimental 
uncertainty. 




° h= 0.002 


* h= 0.003 


+ h = 0.004 


Eq. (4-9) 


Eq. (4-9) 


— - Eq. (4-9) 



Figure 5-1. Ethylene Glycol Temperature Profile 

Flow rate, h=0.002: 77 ml/min 

h=0 . 003 : 153 ml/min 

h= 0.004: 371 ml/min 

Supply Pressure, h=0.002: 40 psi 

h=0 . 003 : 25 psi 

h=0.004: 25 psi 



Figure 5-2 clearly shows that the temperature 
difference at higher RPM for the graphite slurry is almost 
double the temperature difference shown in Figure 5-1 at 
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a h= .002 


* h= .003 


x h= .004 


Eq. (4-19) 


Eq. (4-19) 


Eq. (4-19) 



Figure 5-2. Graphite/Ethylene Glycol Temperature 
Profile 

♦ Flow rate, h=0.002: 46 ml/min 

h=0.003: 114 ml/min 

h=0.004: 136 ml/min 

♦ Supply Pressure, h=0.002: 40 psi 

h=0.003: 30 psi 

h=0.004: 15 psi 



equal RPM and film thickness for the Newtonian fluid, 
ethylene glycol. It is especially prevalent for the 0.002 
inch film thickness. This demonstrates that the presence of 
the graphite powder, thus making the fluid non-Newtonian, 
greatly increases the temperature change. For a rotating 
bearing with microstructures in the lubricant, the increase 
in temperature and subsequent power loss is expected due to 
an increase in the shear stress (15, p. 106) . As discussed 
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in Chapter 4, the Power Law constants are dependent upon 
temperature. The values used for the rotation speed up to 
500 RPM corresponded to temperature readings within the 30°C 
range. When the rotation exceeded 500 RPM, the temperature 
readings were higher and fell within the 40 °C range. 
Therefore, the Power Law constants were adjusted resulting 
in the discontinuity shown in the theoretical curve in 
Figure 5-2 at h=0.002 inches. The uncertainty for the data 
points was too small to plot. The uncertainty for the 
mathematical model was determined as discussed in Chapter 3 
with equation (3-2) . 

Pressure Results 

The purpose of the pressure data was to determine the 
relationship between the pressure profile and the rotation 
speed (RPM) . The inertial and temperature effects on the 
pressure are assumed to be negligible when compared to the 
shear stress. It has already been shown that the addition 
of the graphite particles increases the pressure at lower 
flow rates thus, increasing the load capacity of the 
hydrostatic step bearing (8, p. 33; 14, p. 16; 15, p. 106; 
16, pp. 32-33) . The effect of rotation on the pressure 
ratio versus the radius ratio is shown in Figures 5-3 
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through 5-8. The measured data compared closely with the 
mathematical predictions as developed in Chapter 4 at the 
lower rotation speeds within the uncertainty. 



0.9-- 
0 . 8 " 

0.7- 
0 . 6 - 
ft: 0.5" 

CL 

0.4” 

0.3- * 

0 . 2 -- 

0 . 1 ” : : : : : : : : 

Oi i i i i r i i i i — ^ 

0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 

r/R 



Eq. (4-26) 


* 100 RPM 


♦ 250 RPM 


X 400 RPM 


A 500 RPM 


° 600 RPM 



Figure 5-3. Ethylene Glycol Pressure Profile 
♦ h=0.002; Flow rate=77 ml/min 

The repeatability was within experimental uncertainty. 
The uncertainty for the data points was too small to plot. 
The uncertainty for the mathematical model was determined as 
discussed in Chapter 3 with equation (3-2) . 

The inlet pressure varied slightly during the course of 
the experimentation and therefore, the pressures were 
nondimens ionali zed by the appropriate supply pressure 
reading resulting in a pressure ratio of 1 at the step of 
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r/R 



Eq. (4*26) 


* 100 RPM 


♦ 250 RPM 


X 400 RPM 


A 550 RPM 


a 670 RPM 



Figure 5-4. Ethylene Glycol Pressure Profile 
♦ h=0.003; Flow rate=153 ml/min 




Eq. (4-26) 


* 100 RPM 


+ 250 RPM 


X 400 RPM 


A 550 RPM 


a 670 RPM 



Figure 5-5. Ethylene Glycol Pressure Profile 
♦ h=0.004; Flow rate=371 ml/min 
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r/R 



Eq. (4-28) 


* 100 RPM 


* 250 RPM 


x 400 RPM 


A 500 RPM 


° 600 RPM 



Figure 5-6. Graphite/Ethylene Glycol Pressure 
Profile 



♦ h=0 . 002 ; Flow rate=46 ml/min 




Eq. (4-28) 


* 100 RPM 


+ 250 RPM 


x 400 RPM 


A 550 RPM 


° 670 RPM 



Figure 5-7. Graphite/Ethylene Glycol Pressure 
Profile 

♦ h=0.003; Flow rate=114 ml/min 
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r/R 



Eq. (4-28) 


* 1 00 RPM 


♦ 250 RPM 


X 400 RPM 


A 550 RPM 


o 670 RPM 



Figure 5-8. Graphite/Ethylene Glycol Pressure 
Profile 

♦ h=0.004; Flow rate=136 ml/min 
the bearing for each RPM value. The remaining pressure 
points were nondimens ionali zed by the same supply pressure. 

The 100 RPM pressure ratios in Figures 5-3 and 5-5 
varied slightly from the mathematical model, but are within 
the uncertainty. As shown in Figure 5-3, the increased 
speed resulted in a slight drop in the interior pressure 
readings whereas the two exit pressures did not vary 
appreciably. The probable reason for this decrease in 
pressure is due to the temperature dependency of the fluid 
viscosity. As noted in the previous results, the 
temperature of the bearing increased due to the raised 
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frictional loss at the higher rotations. As temperature 
increases, the viscosity becomes smaller resulting in a 
smaller pressure ratio. This decrease in pressure is also 
indicated in Figures 5-4 and 5-5. Preliminary temperature 
measurements of the interior of the bearing showed higher 
temperatures than those observed at the inlet or outlet 
points of the fluid flow emphasizing the presence of a 
greater temperature affect on the interior pressure reading 
than the outlet readings. The non-Newtonian pressure 
behavior shown in Figure 5-6 also indicates a drop in the 
pressure ratio. This decrease in the pressure ratio is 
larger than shown for an equal film thickness and RPM of the 
ethylene glycol in Figure 5-3. These results are further 
supported by comparing Figures 5-7 and 5-8 to their ethylene 
glycol counterparts, Figures 5-4 and 5-5. As shown in Table 

4- 1, the Power Law constants for the graphite/ethylene 
glycol slurry are also temperature dependent. Therefore, 
the non-Newtonian graphite slurry is also temperature 
dependent. The graphite slurry demonstrated a higher 
pressure drop than the Newtonian fluid, ethylene glycol, 
because of the higher temperatures shown in Figures 5-1 and 

5- 2 and discussed previously. 



41 



The decrease in pressure is significant because it 
translates to a loss of load carrying capability of the 
bearing. The higher load carrying capacity of the 
non-Newtonian lubricant is one of the reasons for their use. 
The pressure is still higher than the Newtonian lubricant, 
but it must be considered because the change in temperature 
experienced at extreme rotations will be substantially 
higher than Newtonian lubricants resulting in a higher 
pressure loss. 



CHAPTER 5 
CONCLUSIONS 



The objective of this research was to determine the 
effects of adding a powder lubricant to a Newtonian carrier 
fluid on the lubrication temperature and pressure profile. 
The effects were shown through result comparisons between 
the pure ethylene glycol and the Graphite/Ethylene Glycol 
slurry. 

A detailed uncertainty analysis of the measuring 
devices was undertaken to determine deviations between 
measured and predicted data. 

Results compared favorably with the mathematical 
models . 

The data presented agreed with the theoretical 
predictions based upon the rheological properties obtained 
from previous papers. These agreements give a high level of 
confidence for the use of the equations derived in this 
paper . 

It was noted that the pressure profile of the graphite 
slurry was affected more by the higher rotations than the 
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pure Newtonian fluid. This affect can be expected due to 
the higher shear rate and subsequent temperature increase. 
Since viscosity of the graphite slurry is temperature 
dependent and the potential applications involve a 
significant increase to the rotation speed, it is 
recommended that further research should be conducted to 
study the interior temperature gains and their affect on the 
pressure profile of the bearing. Measuring the temperatures 
of the interior of the bearing in conjunction with the Power 
Law constants' curve fit as functions of temperature could 
result in a complete understanding of the pressure profile 
of a hydrostatic step bearing subjected to a rotation. 

Preliminary testing showed that the temperature inside 
the bearing is higher than those noted at the inlet or 
outlet of the fluid flow. Therefore, any attempt to 
mathematically model the temperature dependency of the 
viscosity or the Power Law constants would have been 
inaccurate . 

Another factor that could be present at significantly 
higher rotation speeds than those tested in this research is 
possible inertial effects of the graphite powder in the 
slurry mixture. That effect was not apparent in this 
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research since an increase in the pressure profile would be 
expected for such an occurrence. 

The advantages of the non-Newtonian lubricant are their 
increased viscosity and subsequent higher load carrying 
capacity and their ability to operate at higher 
temperatures. However, when the bearing is operated at high 
rotation speeds, the higher shear rate creates an increased 
temperature which could have a weakening affect on their 
pressure profile. This affect would limit their load 
carrying capacity, although it still exceeds the capacity of 
the Newtonian lubricant. This thesis has shown the 
fundamental lubrication principles of a graphite slurry 
non-Newtonian fluid that was pressure induced into a 
rotational hydrostatic step bearing. Further investigations 
into higher rotation speeds and interior temperature gain 
are recommended to more fully understand the non-Newtonian 
behavior of the graphite/ethylene glycol lubricant. 
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APPENDIX SAMPLE UNCERTAINTY CALCULATION 



Determining the uncertainty in the film thickness, 
equation (3-9) , was accomplished through the use of equation 
(3-10) which is expanded as follows: 



CO h = 



dh 



\ 2 



CO 



+ 1 #-C0, 



\ 2 



9yA W »J ^ J 



+ 



dh „ 

dy c ® yc 



21 



1 

2 



Substituting the values 'of the partial derivatives of 
equation (3-9) and the uncertainties of each LVDT results: 

co * = [(-. 0022 (±. 04)) 2 + (-. 0025 (±. 07)) 2 + (~. 0033 (±. 05 )) 2 ] ’ 

0 )* = ± 0.00026 = ± 0.0003 

This is the uncertainty stated on page 16 and 
demonstrates the procedure followed for all the uncertainty 
calculations . 
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